arXiv:1505.00850v2 [cs.IT] 29 Jul 2015 


Full-Duplex Relaying in MIMO-OFDM Frequency-Selective 
Channels with Optimal Adaptive Filtering 

Joao S. Lemos*, Francisco A. Monteiro^^, Ivo Sousa*, and Antonio Rodrigues* 

*Instituto de Telecomunicajoes, and Instituto Superior Tecnico, Universidade de Lisboa, Portugal 
l^lnstituto de Telecomunicafoes, and ISCTE - Instituto Universitario de Lisboa, Portugal 
email: joao.sande.lemos@tecnico.ulisboa.pt , francisco.monteiro@lx.it.pt 


Abstract —In-band full-dnplex transmission allows a relay station to 
theoretically double its spectral efficiency by simultaneously receiving 
and transmitting in the same freqnency band, when compared to 
the traditional half-duplex or out-of-band full-dnplex connterpart. 
Conseqnently, the induced self-interference suffered by the relay 
may reach considerable power levels, which decreases the signal-to- 
interference-plus-noise ratio (SINK) in a decode-and-forward (DF) 
relay, leading to a degradation of the relay performance. This paper 
presents a techniqne to cope with the problem of self-interference 
in broadband mnitiple-input multiple-ontpnt (MIMO) relays. The 
proposed method uses a time-domain cancellation in a DF relay, 
where a replica of the interfering signal is created with the help of a 
recursive least squares (RLS) algorithm that estimates the interference 
frequency-selective channel. Its convergence mean time is shown to 
be negligible by simulation results, when compared to the length 
of a typical orthogonal-frequency division multiplexing (OFDM) 
sequences. Moreover, the bit-error-rate (BER) and the SINR in a 
OFDM transmission are evaluated, confirming that the proposed method 
extends significantly the range of self-interference power to which the 
relay is resilient to, when compared with other mitigation schemes. 

Index terms - Decode-and-forward relay, in-band full-duplex, broad¬ 
band MIMO, freqnency-selective channel, adaptive filtering, recursive 
least squares, self-interference cancellation. 

1. Introduction 

A relay station is a key element in a wireless multi-hop network, 
which is believed to incorporate future communication systems, since 
it may provide a wider coverage, a higher data rate and a lower 
transmit peak power. Therefore, a considerable amount of research 
has been conducted in this topic, where the integration of multiple- 
input multiple-output (MIMO) techniques at the relay appear as a 
natural solution to avoid the key-hole effect (T). 

In-hand full-duplex operation is a novel technique, which has 
recently gained attention in the wireless communication field 
Its main advantage is the possibility to double the spectral efficiency 
of the relay station. This is achieved hy employing the same time 
and frequency resources for receiving and transmitting, i.e., the relay 
operation of receive and forward information is performed within the 
same frequency hand at the same time. However, the clear limitation 
of full-duplex operation is the arising loopback self-interference due 
to the leakage of the relay outgoing signal to the relay receiver 
side, mainly enhanced by the high power unbalance between the 
desired and the self-interference signal, causing inadmissible levels 
of interference and deteriorating the relay performance 0 Thus, the 
self-interference is typically attenuated by a proper system physical 
design and by subtracting a delayed version of the own signal at an 
analog stage (4). Then, a digital baseband processing stage is often 
introduced 0,0. m, so that a reliable communication link may 
be established. 

The relaying protocol often employed in full-duplex systems is 
characterized by regenerating the message from the original source, 
known as decode-and-forward (DF) operation. When compared to 


amplify-and-forward (AF), that blindly repeats the ongoing message, 
DF additional complexity grants an improvement in terms of perfor¬ 
mance 0 Among these digital DF techniques stand two methods, 
time-domain cancellation, subtracting a baseband version of the self¬ 
interference signal, and interference suppression, which exploits the 
degrees of freedom available in MIMO channels 0, 0, Oni. The 
first method chiefly suffers from self-interference erroneous channel 
estimation, while the latter is affected by the distortion from spatial 
shaping of the relay transmitted signal. Besides the aforementioned, 
both performances are deteriorated by limited dynamic range and 
transmission impairments that introduce additional noise sources 
CD. Those effects can drastically harm the operation of full-duplex 
relaying, boosting the search for better mitigation schemes. 

Adaptive cancellation was firstly proposed in im, where the 
authors make use of a gradient-descent algorithm to estimate the 
self-interference and achieve a 28.6 dB of additional interference 
attenuation. This paper proposes for the first time, to the best of our 
knowledge, the use of recursive least squares (RLS) filtering to further 
improve that resilience to self-interference, which converges to the 
optimal mean-square error (MSB) estimator in negligible time. The 
well-known RLS algorithm 1731 is adapted to the MIMO full-duplex 
relay. Moreover, it is a suitable technique to cope with frequency- 
selective effects present in broadband MIMO transmissions, by 
exploiting orthogonal-frequency division multiplexing (OFDM) trans¬ 
missions d. Moreover, this scheme acts as an add-on block to the 
half-duplex relay, i.e., it does not change the relay protocol and can 
be added whenever the full-duplex mode is enabled. 

11. System Model 

Consider a single-hop wireless communication through a MIMO 
relay station, which operates in a decode-and-forward full-duplex 
mode. The relay serves a source with a total of Ns antennas that 
intends to establish a communication link to a destination with Nd 
antennas. The relay is equipped with Mr receiving antennas and 
Mt transmitting antennas. Figure[T]depicts the considered broadband 
uncoded MIMO transmission through a full-duplex relay with self¬ 
interference cancellation. 

At time instant n, a source signal x(n) G is transmitted via a 
relay station to the destination. The relay station receives a version of 
that signal, q(n) G after passing through a frequency-selective 

channel and being mixed with the self-interference and additive noise 
components. At the same time, the relay regenerates and retransmits a 
signal, t(n) G C^®’, which reaches the destination as y(n) G 
also experiencing channel and noise effects. Equation Q shows the 
expressions for both q(n) and y(n) 

q(n) =HsR(2)x(n) -f HLi(z)t(n) -f nR(n), 

y(n) =HRE)(2)t(n) -f nE)(n), 



Figure 1: Model of a decode-and-forward in-band full-duplex relay station with adaptive self-interference cancellation. 


where Hsr( 2) G is the Lsr order channel from the 

source to the relay, Hrd( 2 :) £ Lro order 

channel from the relay to the destination, and Hli( 2) G 
is the Lli order self-interference channel. Both frequency- 
selective channels matrices are represented by their a-transform, 
ii{z) =, where a represents the discrete de¬ 
lay operator and L the filter order; e.g., HsR(z)x(n) = 

vectors nR(n) and nD(n) are the 
additive noise components at the input of the relay station and 
destination, respectively. 

The relay DF protocol /r(-) is independent of the self-interference 
cancellation architecture, and regenerates a delayed estimation of the 
received data streams, t(n) = fr{e{n — d), e{n — d — 1), • • • , e(n — 
d — D)), where d stands for the necessary processing delay and D + 1 
is the length of the employed time-window. Further, the processing 
delay is considered strictly positive and sufficiently long (5), such that 
x(n —fc) and t{n — l) become uncorrelated, i.e., E{x(n —fc)t^(n— 
Z)} = 0, for all fc = 0, ■ • ■ , Lsr and all Z = 0, ■ • ■ , Lli- 

The practical implementation impairments in the relay, arising from 
imperfect AD/DA conversion, limited dynamic range, quantization 
errors, power amplifier non-linearities, etc., are jointly modeled as 
additive interference noise CD, giving t(n) = t(n) + £t{n), where 
t(n) is the known baseband signal at the relay. The error £t(n), takes 
into consideration the relay transmit power weighted by factor S, and 
is characterized by 

£:t(n) ~CAf(0,(5-E{t"(n)t(n)})I). (2) 


Moreover, the relay may use an estimation of the self-interference 
channel matrix, Hli(z), when necessary. However, this estimation 
is associated with an additive noise term, EHLi[fc](n), which models 
the practical non-ideal channel estimation in such relay systems 

HLilk] = ilLi[k]+£H^,lk]{n), for Zc = 0, ■ • ■ , Lu. (3) 
III. Mitigation of Self-Interference 

In order to properly mitigate the self-interference present in the 
relay pre-decoding signal e(n), we apply a feedback canceler filter in 
the time domain. The idea is to estimate the self-interference channel 
effect on the transmitted signal, given by f(n) = HLi(z)t(n), 
employing a finite impulse response (FIR) filter A(z) G 
with order La, that simulates the self-interference chain effect, from 
the relay transmit side to the receive one. Therefore, we have 

e(n) = q(n) + z(n) 

= HsR(2:)x(n) -f HLi(z)t(n) + nR(n) + A{z)t{n). 


A. Time-Domain Cancellation 

Time-domain cancellation (TDC) is the trivial solution to the 
problem presented above, as it uses the self-interference channel 
estimation to perform the cancellation 15), as 

Atdc{z) = -Hli(2), (5) 

where the filter order should be La = Lli- However, this procedure 
is severely harmed by the above mentioned system impairments, as 
we cannot remove the following residual interference 

f(n) + z(n) = HLi( 2 :)t(n) + ATDc{z)l{n) 

= £iiLi{z)(t{n) -\-£tin)'^ + ilLi{z)£t{n), 

being this technique performance dependent on the power of the 
above residual terms. 

B- Recursive Least Squares Cancellation 

The proposed cancellation scheme is based only on the transmitted 
signal t(n) and observed signal q(n), easily accessed by the relay, 
to estimate the self-interference effect. We start by considering the 
mean square error of the self-interference channel estimation effect at 
time instant n, MSE{n, A{z)) = E{(f(n) — A( 2 )t(n))^(f(n) — 
A(«)t(n))}, for the generalized MIMO case, and defining the 
recursive least squares approximation of it as 

MSERLs{n,A{z)) = 

n /x\ 

= - A( 2 )t(fc))^(f(fc) - A{z)i{k)), 

k = l 

where 0 < A < 1 is the forgetting factor, which regulates the de¬ 
pendency on previous observations. The RLS algorithm is derived by 
exploiting the second derivative of the MSErls, as a Newton-type 
algorithm CD . Defining the concatenation matrix of the filter parame¬ 
ters A*(n) = [A[0](n),--- ,A[LA](n)]^ € £^M-r{LA+i)xMfi^ 
RLS correlation matrices as 

n 

fc = l 

n 

= ^A""''Rj(Zc),G c(^A-hl)MTX(LA-hl)MT^ 
k=l 

where Rt(^) = t(n)t^(n), with t(n) = vec{T(n)} = 

vec{ [t(n), • • • , t(n — La)] } G £(^a-i-i)aItxi ^ column vector 
containing the elements of T(n), and (g) represents the Rronecker 


























product, it is possible to write the Newton algorithm time update 
rule by taking the first and second derivatives of equation l|^, as 

A*,„ = - I A*(n - 1))] 

^M^RLs{n,A4n-l)) 

(7) 

where /r is the Newton-type algorithm step size. 

However, equation 0 still depends on f(n), which is in fact 
not available at the relay station. Further, the algorithm requires at 
each iteration the inversion of matrix costing 0{{{La -f 

IV. Algorithm Analysis 

In this section we provide an efficient update rule for the proposed 
algorithm, using the available relay input signal, q(n). Also, the 
asymptotic error covariance matrix is evaluated. 

A. Update Rule 

Firstly, we notice that the correlation between q(n) and t(n), 
assuming weakly stationary variables, is given by 

R-q,t(fe) = E{q(n)t"(n-fc)} 

= E{(HsR( 2 ;)x(n) + f(n) + nR(n))t"(n - fc)} (8) 
= E{f(n)t^(n — fc)}}, 

where E{x(n)t^(n — A:)} = 0, for a sufficiently large processing 
delay, and also E{nR(n)t^(n — k)} = 0. This result allows to 
observe q(n) instead of f(n), since the correlation matrix in 0 is 
~ t(^)’ therefore, not changing the filter convergence 

properties. Nevertheless, the equivalent noise present in the observed 
signal q(n) is now composed by the source and additive noise vec¬ 
tors, which will affect the filter performance, mainly its convergence 
time. 

Resorting to the Woodbury matrix identity ini, by defining 
P(n) = it is possible to show that 

P(n + 1) =[AE?E^^(n - 1) + R,-(n)] 

= i[p-i(n) + iRE(n)]-' 

= i[P~^(n)+ t(n)it^(n)]"^ 

- P(^)R-t(w)P(n) \ 

AV '' '' A-f t»(n)P(n)t(n)/’ 

which does not require the matrix inversion in 0, reducing the 
update rule complexity order from cubic to quadratic. 

Thus, we can finally present the algorithm 3-step update rule as 

kfni = _P(n)t_(n) _ 

A-f t^(n)P(n)t(n) ’ 

L‘A 

A*,„ = A*,„_i -b/rk(n)(q(n) -^A[/](n- l)t(n-/))'^, 

1=0 

P(n-l-l) = i ^P(n-b 1) — k(n)t^(n)P(n)^, 

( 10 ) 


B. Asymptotic Error Covariance 


In this subsection, we prove that the algorithm converges to the 
optimal value by evaluating the asymptotic error covariance matrix of 
the iterative expressions in l |10| l. For that end, we assume the case of 
no forgetting factor (A = 1), which guarantees the filter convergence 
1151 . Considering the Newton update rule in l[7} and the observation 



Thus, as n grows large, the RLS covariance matrices will converge 
to the true covariance matrices, ®(n) —> Sjj = E{t(n)t^(n)} 
and > Sj ^ = E{q^(n) (g) t(fc)}, by the law of large 

numbers O- We may then write the error matrix at instant n as in 



— ■^-k,Opt. 


= [S?^^®(n)] (n)-EriE, 


■'t.t t,q 


= [S?t-®(«)] 'y]q"(fe)®t(fc)-A*,opt. (11) 

A: = l 



(g)t(fc), 

k = l 


where the last step comes from the above approximation, some 
algebraic manipulation, and defining f{k) — q_{k) — Aopt.{z)t{k) 
as the estimation error vector. Finally, we reach the covariance error 
parameter given by 

Q(n) = E{A*,„A^^„} 

1 n n ^ 

k=l 1=1 

n n 

- Sj-i^E{y]y]t(fe)F(fc)f(/)t^(Z)}Eri 

k=l1=1 
n 71 

k=l1=1 

( 12 ) 


where it was assumed that f(n) and t(n) are mutually independent. 
Thus, the error power of the RLS estimation for A = 1 converges to 
zero as the number of iterations n tend to infinity, which provides an 
optimal estimation of the self-interference matrix. 


V. Numerical Results 

In this section, we study the proposed RLS filter performance, 
in terms of its convergence time, SINR and bit-error-rate (BER), 
evaluated at the relay station. All results are obtained via Monte- 
Carlo simulation of 2000 OFDM transmitted symbols, for a typical 
relay station scenario. 


A. System Parameters 

We consider an uncoded MIMO-OFDM system, where Ns = 
2 data streams transmit a 16-QAM modulated OFDM stream of 
symbols with a block of Nsub subcarriers and a cyclic prefix of 
Acp = Lsr, to a destination also with No ~ 2 antennas. The relay 
is considered symmetric and composed by Mr = 3 receive antennas 
and Mt = 3 transmit antennas, with a proper rate adaptation. Further, 
the channels are assumed to be of order one, i.e., Lsr = Lrd = 
Lri = 1 (2 taps channels), similar to what is considered in 1121 . The 
channels Hsr( 2) and Hre)( 2 ) are drawn from complex Gaussian 
distributions and taken from CN{Q, 1), while each self-interference 
matrix channel tap has distribution HLi[fc] ~ CA/'(0, for 


where vector k(n) is the update direction of the filter 1131 . 









Figure 3: Impact on the SINK of the self-interference power for the 
RLS, TDC and NI methods. 


Figure 2: Histogram of the proposed RLS algorithm convergence time 
for the self-interference channel estimation. 


fe = 0, • • • , Lli, where accounts for the residual power of the 
self-interference channel, after propagation and analog-circuit first 
stage of mitigation. The filter order is set to be La — Lli = 1, 
so that we can achieve a perfect estimation of the self-interference 
channel. The transmitted source and relay signals have normalized 
power, i.e., before OFDM modulation we have E{x^(n)x(n)} = 1 
and E{t^(n)t(n)} = 1. Additionally, the thermal noise present 
at the relay, nR(n), is considered to be additive white Gaussian 
with distribution nR(n) ~ CA/'(0, where its power is set 

to = —15 dB. Finally, the filter update rule is as defined in 
equation Co} with no forgetting factor (A = 1) and with step size 
fi=l. 

B. Convergence Time 

Let us firstly analyse the convergence time of the RLS algorithm, 
i.e., the number of iterations required from an initial starting point, 
= 0 and P(l) = I, to a certain point where the filter parameters 
satisfy the following error metric (EM) C3 

(13) 

II -H-li,* II p 

where Hli.* is the matrix concatenation of the self-interference 
channel taps, Hli,* = [Hli[0], • ■ • ,Hli[La]]^. 

To that end, we consider an OFDM transmission with Naub = 
8192 subcarriers and a self-interference channel with power (t|,/ = 0 
dB. The implementation impairments are taken into consideration by 
setting (5 = 10“®, as in 0. Figure shows the histogram of the 
self-interference channel estimation RLS algorithm convergence time 
distribution, considering EM < —30 dB, after 20000 realizations. 
Further, a superimposed log-normal distribution indicates that the 
mean value of the algorithm convergence time is 1007 samples. This 
number is in fact less than the OFDM symbol length considered for 
simulation, and negligible among the number of transmitted OFDM 
sequences in real implementations. 


C. SINK and BER 


We now evaluate the SINR and the BER for the parameters 
considered above and for different values of interference power 
o-Lj. We compare the proposed cancellation technique with the case 
where there is no filtering. A* = 0, known as natural isolation 
(NI), and TDC with a channel estimation error of —20 d B, i.e. 
£jri^j[k]{n) ~ CAf{0, aaLfI) with a = 10“^, as in section 


III-A 



Figure 4: BER comparison of the proposed RLS algorithm, TDC and 
NI, for a 16-QAM OFDM zero-forcing detector. 


The BER gain against an equivalent half-duplex (HD) system is also 
assessed. Eigurej^shows the SINR curves, defined as Px/ (Pi+Pnji), 
where i(n) = f(r!,) + z(n) is the remaining interference component 
after cancellation, for the three considered cases, when the RLS 
algorithm has converged. There, we can observe the large gain 
provided by the RLS cancellation, which is not affected by erroneous 
channel estimation as the TDC case. The provided SINR gain in 
Eigure consequently translates into a BER curve gain, as shown 
in Eigure By applying the proposed RLS filter, it is possible to 
obtain a gain of around 15 dB of self-interference resilience when 
compared to the case of TDC. Altogether, approximately 35 dB are 
obtained when taking the NI performance as the base level, which 
surpasses the 28.6 dB margin obtained in OH. 

VI. Conclusions 

It was shown that the multi-dimensional RLS filter is capable of 
efficiently estimate the self-interference in an in-band full-duplex 
decode-and-forward broadband MIMO relay. The paper first derived 
the MIMO RLS filter in the context of full-duplex relaying and then 
assesses its performance in comparison to other traditional mitigation 
techniques. A computationally efficient filtering update rule was 
defined for this application and it was latter shown that it assures 
convergence to the optimal MSE estimation. The convergence time of 
the algorithm was numerically evaluated and shown to be negligible 
when compared to a practical wireless OEDM transmission block. 
Einally, the simulation of a full-duplex relay system allowed us to 
evaluate the interference resilience of the proposed technique, in 
terms of SINR and BER. In the case of both these metrics, one can 
found a very significant gain when compared to previously proposed 
signal processing-domain mitigation techniques. 
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